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STATEMENT  OF  PROBLEM  STUDIED 


The  research  project  described  here  was  concerned  with  the  develop¬ 
ment  of  multiaxial  constitutive  equations  for  low-to-medium  density  snow. 
The  constitutive  equations  were  to  have  the  following  characteristics: 

(1)  describe  the  material  behavior  in  terms  of  the  microstructural 
properties  of  the  material 

(2)  represent  the  material  for  intermediate  and  high  strain  rates, 
large  strains,  and  multiaxial  states  of  stress. 

The  constitutive  formulation  would  have  to  satisfy  the  requirements  of 
material  frame  indifference  and  have  a  thermodynamic  basis  in  that  it 
would  be  consistent  with  the  second  law  of  thermodynamics. 

The  first  criterion  listed  above  would  require  a  major  development 
in  itself,  since  measuring  important  microstructural  properties  is  a 
difficult  problem.  A  method  for  (i)  determining  which  microstructural 
properties  (such  as  grain  size,  pore  size,  bond  length,  bond  radius, 
density,  bonds/grain,  specific  surface  area,  etc.)  are  relevant  to 
deformation,  (ii)  quantitatively  measuring  these  parameters,  and  (iii) 
determining  how  they  vary  with  macroscopically  observable  variables  such 
as  strain  and  strain  rate. 

In  order  to  accomplish  this,  a  complete  stereological  theory  had  to 
be  formulated.  Current  stereological  formulations,  to  the  author's 
knowledge,  were  not  adequate  for  accurately  determining  these 
microstructural  parameters  for  porous  materials  with  irregular  grain 
size,  shape,  and  distribution. 

Once  the  stereological  analysis  was  developed,  a  constitutive  theory 
then  had  to  be  considered.  During  the  past  thirty  to  forty  years,  a 
number  of  formulations  have  been  used  to  describe  the  behavior  of  snow. 


These  formulations  have  included  linear  elastic,  linear  viscoelastic, 
nonlinear  viscoelastic,  Newtonian  fluid,  nonNewtonian  fluid,  and  elas¬ 
tic-plastic  models.  None  of  these  models  can  be  considered  to  be  truly 
correct,  since  snow  is  neither  elastic,  viscoelastic,  fluid  or  plastic  in 
behavior.  However  any  of  the  above  formulations  can  give  a  good 
approximation  of  the  material  behavior  if  certain  conditions  (usually 
quite  restrictive)  on  the  deformations  are  met. 

What  was  to  be  attempted  in  this  project  was  the  development  of  a 
constitutive  equation  which  would  enable  one  to  determine  alterations  in 
the  microstructure  of  the  material  due  to  loading.  This  constitutive  law 
would  also  be  more  general  than  preceding  formulations  and  hence  could  be 
applied  to  a  broader  class  of  problems.  In  addition  it  was  anticipated 
that  this  formulation  would  add  to  current  understanding  of  the  important 
deformation  mechanisms.  Once  the  general  formulation  was  finished,  it 
was  also  intended  to  consider  simplified  forms  which  could  be  used  for 
special  situations. 

Finally,  some  testing  was  to  be  done  in  order  to  evaluate  the  above 
results.  A  stress  wave  test  program  was  to  be  instituted  to  observe  the 
alteration  of  the  material  microstructure  due  to  shockwaves.  Also  the 
results  would  be  compared  to  those  predicted  by  the  constitutive  theory. 
SUMMARY  OF  RESULTS 

The  research  project  was  divided  into  three  main  parts,  based  on  the 
work  effort  needed  to  achieve  the  project  goals.  These  are:  (a) 
development  of  a  stereological  theory  for  evaluating  the  microstructure 
of  snow,  (b)  formulation  of  a  constitutive  equation,  and  (c)  a  stress 


wave  test  program. 


The  first  part  of  the  project  required  a  major  development  effort, 
since  current  stereological  theory  was  not  developed  to  the  point  where 
it  could  be  used  to  make  a  detailed  microstructural  analysis  of  a  materi¬ 
al  with  highly  irregular  grain  structure.  The  ability  to  measure  such 
properties  as  mean  grain  size,  pore  size,  density,  bond  diameter,  neck 
length  and  the  number  of  bonds  per  grain  was  mandatory.  These  parameters 
were  to  be  used  in  the  formulation  of  a  constitutive  theory  for  snow, 
since  they  were  considered  to  be  the  important  microstructural  properties 
which  determine  the  mechanical  properties  of  snow.  These  parameters  were 
to  be  used  in  constructing  the  constitutive  theory  for  snow. 

A  stereological  formulation  was  developed.  This  theory  allows  one 
to  calculate,  with  the  use  of  a  surface  section  of  snow,  all  of  the  above 
mentioned  microstructural  parameters.  These  calculations  can  be  made 
irregardless  of  the  grain  size  and  shape  distributions.  To  the  author's 
knowledge,  this  represents  a  new  development.  It  also  enables  one  to 
avoid  using  thin  sections  to  analyze  the  microstructure  of  snow.  Thin 
section  techniques  entail  a  more  tedious  procedure  to  acquire  the  same 
information. 

The  second  portion  of  the  project  was  the  development  of  a 
multiaxial  constitutive  theory  for  large,  high  rate  deformation  'f  snow. 
This  theory  models  the  material  as  a  visco-plastic  material.  The 
Lagrangian  strain  rate  is  determined  to  be  a  function  of  the  stress, 
stress  rate,  and  a  state  vector  which  represents  the  microstructural 
properties  of  the  material.  The  history  dependence  of  the  material  is 
defined  in  terms  of  evolution  relations  which  determine  the  rate  of 
change  of  the  state  variables  in  terms  of  the  current  state  variables  and 
the  stress  state.  The  state  vector  has  as  components  the  density,  grain 
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size,  pore  size,  bond  diameter,  neck,  length,  number  of  bonds  per  grain, 
and  the  surface  area  per  unit  volume  of  the  material . 

The  resulting  formulation  goes  a  long  way  toward  defining  much  more 
accurately  then  previously  done  the  relative  importance  of  the  various 
deformation  mechanisms  which  contribute  to  the  overall  deformation  of  a 
porous  material  such  as  snow.  However,  it  is  a  complicated  theory 
requiring  a  computer  program  to  use  the  formulation.  The  next  step  in 
this  development  would  be  a  program  to  specialize  this  formulation  to 
special,  restricted  classes  of  deformation  in  which  only  one  or  two  of 
the  various  mechanisms  are  predominant.  In  this  way  valid  simplified 
equations  can  be  generated. 

The  final  stage  of  the  project  was  a  stress  wave  testing  program. 
This  was  carried  out  during  the  1984-85  winter.  This  program  was  run  in 
order  to  provide  means  of  evaluating  the  constitutive  equation  and  to 
learn  more  about  stress  wave  propagation  in  snow.  The  microstruccure  of 
the  material  was  evaluated  both  before  and  after  the  stress  wave  shots. 
Pressure  profiles,  attenuation  rates  and  air  wave  profiles  were  all 
measured  in  the  tests.  Unfortunately,  only  a  limited  number  of  tests 
were  completed  before  two  of  the  three  pressure  transducers  were  de¬ 
stroyed.  However,  the  acquired  data  does  appear  to  be  sufficient  to 
evaluate  the  constitutive  formulation.  This  evaluation  process  is  still 
in  progress  and  should  be  finished  sometime  in  October. 

Finally,  one  project  was  added  to  the  objectives  of  the  proposal. 
This  involved  the  development  of  a  microprocessor-based  penetrometer  to 
quickly  and  accurately  obtain  a  permanent  record  of  the  strength  profile 
of  snowcover.  This  instrument  was  of  considerable  use  in  finding  uniform 


snow  layers  in  which  to  conduct  the  stress  wave  experiments.  However, 


the  instrument  has  also  sparked  much  interest  among  people  who  are 
involved  in  avalanche  hazard  forecasting,  since  it  provides,  in  five 
minutes,  more  information  than  the  Ram  penetrometer  provides  in  an  hour. 
The  Ram  penetrometer  is  currently  the  standard  instrument  used  in  ava¬ 
lanche  hazard  evaluation. 
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